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Abstract 

We present an improved measurement of the standard model CP violation parameter sin2^i 
(also known as sin 2(3) based on a sample of 85 x 10 6 BB pairs collected at the T(45) resonance 
with the Belle detector at the KEKB asymmetric-energy e + e~ collider. One neutral B meson 
is reconstructed in a J/ipKg, ip(2S)Kg, Xc\K% VcKs^ J/^K* , or J/^Kl CP-eigenstate decay 
channel and the flavor of the accompanying B meson is identified from its decay products. From 
the asymmetry in the distribution of the time interval between the two B meson decay points, we 
obtain sin2^i = 0.719 ± 0.074(stat) ± 0.035(syst). 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw 
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In the standard model (SM), CP violation arises from an irreducible complex phase in 
the weak interaction quark-mixing matrix [Cabibbo-Kobayashi-Maskawa (CKM) matrix] [[!]. 
In particular, the SM predicts a CP- violating asymmetry in the time-dependent rates for 
B° and B° decays to a common CP eigenstate fcp, where the transition is dominated by 
the b — > ccs process, with negligible corrections from strong interactions 0: 

A i+\ — — » fcp) - r(.B — > fcp) . . , . m 

A{t) = rffl ^ f CP ) + r ( B° ^ / CP ) = " f ' sm 2 ^ 8m(A " a) ' (1) 

where 1(5°, ZJ° -> / CP ) is the rate for 5° or W to / CP at a proper time t after production, 
£/ is the CP eigenvalue of fcp, Am^ is the mass difference between the two B° mass 
eigenstates, and 0i is one of the three interior angles of the CKM unitarity triangle, defined 
as 0i = ii — a.Tg(y t * b Vtd/V* b V c d). Non-zero values for sin20i have been reported by the Belle 
and BaBar groups [|. 

Belle's published measurement of sin 20! is based on a 29.1 fb _1 data sample containing 
31.3 x 10 6 BB pairs produced at the T(AS) resonance. In this paper, we report an im- 
proved measurement that uses 85 x 10 6 BB pairs (78 fb _1 ). Two important changes exist 
in the analysis with respect to the published result ||; we apply a new track reconstruction 
algorithm that provides better performance and a new proper-time interval resolution func- 
tion || that reduces systematic uncertainties in sin20x. The data were collected with the 
Belle detector || at the KEKB asymmetric collider Jf], which collides 8.0 GeV e~ on 3.5 GeV 
e + at a small (±11 mrad) crossing angle. We use events where one of the B mesons decays 
to fcp at time tcp, and the other decays to a self-tagging state / tag , which distinguishes 
B from B°, at time itag- The CP violation manifests itself as an asymmetry A(At), where 
At is the proper time interval between the two decays: At = tcp — t tSug . At KEKB, the 
T(4S f ) resonance is produced with a boost of f3~f = 0.425 nearly along the z axis defined as 
anti-parallel to the positron beam direction, and At can be determined as At ~ Az/(f3j)c, 
where Az is the z distance between the fcp and / tag decay vertices, Az = zqp — z tSug - The 
average value of Az is approximately 200 /im. 

The Belle detector |j| is a large-solid-angle spectrometer that includes a silicon vertex 
detector (SVD), a central drift chamber (CDC), an array of aerogel threshold Cerenkov coun- 
ters (ACC), time-of-flight (TOF) scintillation counters, and an electromagnetic calorimeter 
comprised of CsI(Tl) crystals (ECL) located inside a superconducting solenoid coil that pro- 
vides a 1.5 T magnetic field. An iron flux-return located outside of the coil is instrumented 
to detect K\ mesons and to identify muons (KLM). 

We reconstruct B° decays to the following CP eigenstates M: J/ipKg, ip(2S)Kg, XciK®, 
rj c Kg for £/ = —1 and J/ipK® for £/ = +1. We also use B° — > J/tpK*° decays where 
K*° —>■ Kg7T°. Here the final state is a mixture of even and odd CP, depending on the 
relative orbital angular momentum of the J/ip and K*°. We find that the final state is 
primarily £ f = +1; the = -1 fraction is 0.19 ± 0.02(stat) ± 0.03(syst) @. J/ip and 
ip(2S) mesons are reconstructed via their decays to (£ = /i, e). The ip(2S) is also 

reconstructed via J/%p7i + Ti~, and the Xd v i a J/ipl- The rj c is detected in the KgK^ii + , 
K + K~ir° J and pp modes. For the J/ipKg mode, we use Kg — > tt + tt~ and 7r°7r° decays; 
for other modes we only use Kg — > 7t + 7t~. For reconstructed B — > fcp candidates other 
than J/ipK^, we identify B decays using the energy difference AE = E™ s — E™^ m and 
the beam-energy constrained mass Mt, c = ^/(-^bcam) 2 — (p C B iS ) 2 i where E^ m is the beam 
energy in the center-of-mass system (cms) of the T(45) resonance, and E c ^ s and p c ™ are 
the cms energy and momentum of the reconstructed B candidate, respectively. Figure [I] 
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FIG. 1: The beam-energy constrained mass distribution for all decay modes other than J/t/jK^ 
(left). The p c ™ distribution for B° — > J/ipK® candidates with the results of the fit (right). 



TABLE I: The numbers of reconstructed B — > fcp candidates before flavor tagging and vertex 
reconstruction, iV rec , the numbers of events used for the sin2<^i determination, A cv , and the esti- 
mated signal purity in the signal region for each fcp mode. 



Mode 





A rcc 


A ev 


Purity 


J/^(£ + £-)K s (tt+tt-) 


-1 


1285 


1116 


0.976 ± 0.001 


J/ij(£ + e-)K° s (-K TT°) 


-1 


188 


162 


0.82 ±0.02 


ifj(2S)(£ + £-)K° s (ir + iT-) 


-1 


91 


76 


0.96 ±0.01 


ifj(2S)(J/ijiT + -K-)K®(iT + ir-) 


-1 


112 


96 


0.91 ±0.01 


Xci(J/^)K° s (7r+7r-) 


-1 


77 


67 


0.96 ±0.01 


ri c {KlK-TT + )K° s {iT + iT-) 


-1 


72 


63 


0.65 ±0.04 


rj c (K + K-n°)K%(Tr+Ti -) 


-1 


49 


44 


0.72 ± 0.04 




-1 


21 


15 


0.94 ±0.02 


All with £f = -1 


-1 


1895 


1639 


0.936 ± 0.003 


J/^(£ + £-)K* (K s tt°) 


-1(19%)/+1(81%) 


101 


89 


0.92 ±0.01 


J/^(£ + £-)K° L 


+1 


1330 


1230 


0.63 ± 0.04 


All 




3326 


2958 


0.81 ±0.01 



(left) shows the M^ c distributions for all B° candidates except for B° — ► J/tpK^ that have 
AE values in the signal region. Table | lists the numbers of observed candidates, N Tec . 

Candidate B° — ► J /ipK° L decays are selected by requiring ECL and/or KLM hit patterns 
that are consistent with the presence of a shower induced by a K\ meson. The centroid of 
the shower is required to be within a 45° cone centered on the K\ direction inferred from 
two-body decay kinematics and the measured four-momentum of the J/ip. Figure [3] (right) 
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TABLE II: The event fractions ei, wrong tag fractions wi, and effective tagging efficiencies e s = 
e\(l — 2wi) 2 for each r interval. The errors include both statistical and systematic uncertainties. 
The event fractions are obtained from the J/ipK s simulation. 



I r interval e/ wi e s 

1 0.000 - 0.250 0.398 0.458 ± 0.006 0.003 ± 0.001 

2 0.250 - 0.500 0.146 0.336 ± 0.009 0.016 ± 0.002 

3 0.500 - 0.625 0.104 0.228 ± 0.010 0.031 ± 0.002 

4 0.625 - 0.750 0.122 0.160 ±Slool 0.056 ± 0.003 

5 0.750 - 0.875 0.094 0.112 ±0.009 0.056 ± 0.003 

6 0.875 - 1.000 0.136 0.020 ± 0.006 0.126 + 2*22f 



shows the p c ^ s distribution, calculated with the B° — > J/ipK® two-body decay hypothesis. 
The histograms are the results of a fit to the signal and background distributions. There 
are 1330 entries in total in the 0.20 < p c ™ < 0.45 GeV/c signal region; the fit indicates a 
signal purity of 63%. The reconstruction and selection criteria for all fcp channels used in 
the measurement are described in more detail elsewhere |J. 

Charged leptons, pions, kaons, and A baryons that are not associated with a reconstructed 
CP eigenstate decay are used to identify the 6-flavor of the accompanying B meson. Based 
on the measured properties of these tracks, two parameters, q and r, are assigned to an 
event. The first, q, has the discrete value +1 (—1) when the tag-side B meson is likely 
to be a B° (B°), and the parameter r is an event-by-event Monte-Carlo-determined flavor- 
tagging dilution factor that ranges from r = for no flavor discrimination to r = 1 for 
an unambiguous flavor assignment. It is used only to sort data into six intervals of r, 
according to estimated flavor purity. The wrong-tag probabilities, wi (I = 1,6), that are 
used in the final fit are determined directly from data. Samples of B° decays to exclusively 
reconstructed self-tagging channels are utilized to obtain wi using time-dependent B°-B° 
mixing oscillation: (Aof — N$p) / (^of + A*sf) = (1 — 2wi) cos(AmrfAt), where Aof and Nsf 
are the numbers of opposite and same flavor events. The event fractions and wrong tag 
fractions are summarized in Table |T[ The total effective tagging efficiency is determined 
to be e cS = Y^=i e i(^ ~ 2u>;) 2 = 0.288 ± 0.006, where ei is the event fraction for each r 
interval. The error includes both statistical and systematic uncertainties. Improvements in 



the Monte Carlo simulation [[10] and in the track reconstruction yield e e g that is higher by 
6.7% (relative) than the value in Ref. ||. 

The vertex position for the fcp decay is reconstructed using leptons from J/tp decays or 
charged hadrons from rj c decays, and that for / tag is obtained with well reconstructed tracks 
that are not assigned to fcp- Tracks that are consistent with coming from a K s — > 7r + 7r~ 
decay are not used. Each vertex position is required to be consistent with the interaction 
region profile, determined run-by-run, smeared in the r-0 plane to account for the B meson 
decay length. With these requirements, we are able to determine a vertex even with a 
single track; the fraction of single-track vertices is about 10% for zqp and 22% for z tag . The 
proper-time interval resolution function _R S j g (Ai) is formed by convolving four components: 
the detector resolutions for zap and z tSug , the shift in the z tag vertex position due to secondary 
tracks originating from charmed particle decays, and the kinematic approximation that the B 
mesons are at rest in the cms |J. A small component of broad outliers in the Az distribution, 
caused by mis-reconstruction, is represented by a Gaussian function. We determine twelve 
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At (ps) 

FIG. 2: The At distributions for the events with q£f = +1 (solid points) and q£f = — 1 (open 
points). The results of the global fit with sin20i = 0.719 are shown as solid and dashed curves, 
respectively. 



resolution parameters from fits of data to the neutral and charged B meson lifetimes || and 
obtain an average At resolution of ~ 1.43 ps (rms). The width of the outlier component 
is determined to be 42^ ps; the fractions of the outlier components are (2 ± 1) x 10~ 4 for 
events with both vertices reconstructed with more than one track, and (2.7 ±0.2) x 10 -2 for 
events with at least one single-track vertex. 

After flavor tagging and vertexing, we find 1465 events with q — +1 flavor tags and 
1493 events with q — — 1. Table | lists the numbers of candidates used for the sin 20! 
determination, N ev , and the estimated signal purity in the signal region for each fcp mode. 
Figure |2] shows the observed At distributions for the q£f = +1 (solid points) and q£f — — 1 
(open points) event samples. The asymmetry between the two distributions demonstrates 
the violation of CP symmetry. We determine sin 20i from an unbinned maximum-likelihood 
fit to the observed At distributions. The probability density function (PDF) expected for 
the signal distribution is given by 

e -\At\/r B o 

V sig (At,q,w h Z f ) = — [1 -qt f (l-2w l )sm2(l) 1 sm(Am d At)}, (2) 

At b o 

where we fix the B° lifetime t#o and mass difference at their world average values [|TTJ. Each 
PDF is convolved with the appropriate i? S i g (At) to determine the likelihood value for each 
event function of sin 20i : 

Pi = (Woi) ^[UzV^At^q^i^^R^At-At') 

+ (1 - U g )V hkg {At')R bkg {At - At')] dAt' + / iP i(At), (3) 

where / S ig is the signal fraction calculated as a function of p c ™ for J j ipK^ and of AE and 
Mbc for other modes. Pbkg(At) is the PDF for combinatorial background events, which 
is modeled as a sum of exponential and prompt components. It is convolved with a sum 
of two Gaussians, iiW, which is regarded as a resolution function for the background. To 
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TABLE III: The numbers of candidate events, 
(statistical errors only). 


N ev , and values of sin20i 


for various subsamples 


Sample 


AW 


sin 2<j)i 


T/ih (V+7r _N i 
J /V JJX sy 7r 71 ) 

{cc)K° s except J/^K^tt+tt-) 
T/il)K* (K ir ) 


1 1 1 

523 
1230 


n 7q 4- n 1 n 

U. I O HI u. ±u 

0.67 ±0.17 
0.78 ±0.17 

U.UI —1— VJ.UO 


/t^ = B°(q= +1) 

J tag ^ \H J- J 


1465 
1493 


0.65 ±0.12 
0.77 ± 0.09 


< r < 0.5 
0.5 < r < 0.75 
0.75 < r < 1 


1600 
658 
700 


1.27 ±0.36 
0.62 ±0.15 
0.72 ± 0.09 


data before 2002 
data in 2002 


1587 
1371 


0.78 ±0.10 
0.65 ±0.11 


All 


2958 


0.72 ± 0.07 



account for a small number of events that give large At in both the signal and background, we 
introduce the PDF of the outlier component, P Q i, and its fraction / [. The only free parameter 
in the final fit is sin20i, which is determined by maximizing the likelihood function L = 
Yii Pi, where the product is over all events. The result of the fit is 

sin 20i = 0.719 ± 0.074(stat) ± 0.035(syst). 

The systematic error is dominated by uncertainties in the vertex reconstruction (0.022). 
Other significant contributions come from uncertainties in w\ (0.015), the resolution function 
parameters (0.014), a possible bias in the sin 20! fit (0.011), and the J/ipK^ background 
fraction (0.010). The errors introduced by uncertainties in Am^ and t b o are less than 0.010. 

Several checks on the measurement are performed. Table [TT1] lists the results obtained 
by applying the same analysis to various subsamples. All values are statistically consistent 
with each other. Figures [3](a)-||(c) show the raw asymmetries and the fit results for all 
modes combined, (cc)Kg, and J/ipK®, respectively. A fit to the non-CP eigenstate modes 
B° — > Z)(*)~7r + , D*~p + , J/ipK* (K + n~), and D*~£ + v, where no asymmetry is expected, 
yields 0.005 ± 0.015(stat). Figure |||(d) shows the raw asymmetry for these non-CP control 
samples. 

The signal PDF for a neutral B meson decaying into a CP eigenstate [Eq. @] can be 
expressed in a more general form as 



e -\At\/r B o 

P sig ( At, q,wi)= { 1 + qil 



2wi 



21mA 
|A| 2 ±1 



sin(AmrfAt) 





A 


2 -l 


+■ 


A 


2 ±1 



cos(Am^At) 



(4) 



where A is a complex parameter that depends on both B°-B° mixing and on the amplitudes 
for B° and B° decay to a CP eigenstate. The presence of the cosine term (|A| ^ 1) would 
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FIG. 3: (a) The raw asymmetry for all modes combined. The asymmetry for J/t/jK^ and J/ipK* 
is inverted to account for the opposite CP eigenvalue. The corresponding plots for (b) (cc)Kg, 
(c) J/ipK®, and (d) non-CP control samples are also shown. The curves are the results of the 
unbinned maximum likelihood fit applied separately to the individual data samples. 
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indicate direct CP violation; the value for sin20i reported above is determined with the 
assumption |A| = 1, as |A| is expected to be very close to one in the SM. In order to test 
this assumption, we also performed a fit using the above expression with acp = — £/ImA/|A| 
and |A| as free parameters, keeping everything else the same. We obtain 

|A| = 0.950 ± 0.049(stat) ± 0.025(syst) 

and acp = 0.720 ± 0.074(stat) for all CP modes combined, where the sources of the sys- 
tematic error for |A| are the same as those for sin 20!. This result is consistent with the 
assumption used in our analysis. 
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